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Abstract. We consider a cosmological model with non- 
Gaussian initial perturbations, which in principle could 
be generated in non-standard inflationary scenarios with 
two or more scalar fields. In particular we focus our at- 



tention on the model proposed by Linde and Mukhanov 



(1997 ), where perturbations are quadratic in a Gaussian 
field. These perturbations, if they exist, have to be observ- 
able as a non-Gaussian distribution of the CMB signal on 
the sky. In order to efficiently pick up the non-Gaussian 
signal in CMB maps of degree resolution, one can use 
Minkowski Functionals and peak statistics. Our paper 
contains the theoretical predictions of the properties for 
Minkowski Functionals and distributions of peaks of the 
CMB anisotropy in the model with "squared" Gaussian 
statistics. Likelihood comparison of the four-year COBE 
DMR data to this non-Gaussian model and the standard 
Gaussian model does not select any of them as most likely. 

We also suggest an efficient algorithm for fast simula- 
tion of CMB maps on the whole sky. Using a cylindrical 
partition of the sphere, we rewrite the spherical harmonics 
analysis as a Fourier transform in flat space, which makes 
the problem accessible to numerically advantageous FFT 
methods. 



the first detection of the CMB anisotropy by the COBE 



satellite ( Smoot et al. 1992 ; Bennett et al. 1996 ) several 



groups have also reported observations on angular scales 
of about 1° which roughly corresponds to the horizon size 
at recombination. Obtaining the spectrum of the primor- 
dial anisotropies on smaller scales with high precision will 
allow us to verify several robust predictions of inflation 
and permit to determine the values of important cosmo- 
logical parameters with high accuracy. Future experiments 
( [Bennett et al. 1995| ; [Bersanelli et al. 1996| ) will yield all- 



sky maps of the CMB with sufficiently high resolution and 
sensitivity. 

Apart from estimates of the power spectrum with sub- 
sequent parameter determination, the CMB data also pro- 
vides very important constraints on the nature of the per- 
turbations that led to the formation of large-scale struc- 



ture (Bond and Jaffe 1998). To extract this information 
beyond the power spectrum we need to consider not only 
the amplitudes, but also the phases of the temperature 
field. One possible approach is under consideration in our 
paper. 

Generally speaking, cosmological models based on the 
inflationary paradigm predict adiabatic Gaussian fluctua- 



Key words: Methods: numerical, Methods: statistical, 
cosmic microwave background, Cosmology: observations 

1. Introduction 

Observations of the Cosmic Microwave Background pro- 
vide fundamental information about the primordial inho- 
mogeneities in the Universe. The temperature anisotropies 
on scales above a few degrees preserve information about 
primordial density fluctuations on scales larger than the 
acoustic horizon at the moment of recombination. Since 
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scale invariant one ( 


Mukhanov and Chibisov 1981 


; see 


also Starobinsky 1982; 


Hawking 1982|; Guth and Pi 1985 


Bardeen et al. 1983 


) . The Gaussianity of the density per- 



turbations directly translates into Gaussianity of the CMB 
temperature fluctuations on the sky. 

However, along with standard inflationary models, at 
present there still exist theories which are compatible with 
observations and predict non-Gaussian primordial fluctu- 
ations. Among them are, for instance, non-standard infla- 
tionary scenarios with two or more scalar fields, where 
one could expect significant deviations from Gaussian- 



ity (Linde and Mukhanov 1997; Antoniadis et al. 1997; 
Peebles 1999a,|hj), namely, the perturbations there are 



quadratic in a Gaussian field. Another possibility to accu- 
mulate non-Gaussianity in the CMB signals exists after 
recombination. Even if the fluctuations were Gaussian at 
the surface of the last scattering they may have acquired 
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non-Gaussian contributions due to subsequent weak grav- also find significant deviations from Gaussianity. It was 



itational lensing (Fukushige et al. 1995; Bernardeau 1997) 
and due to various foregrounds like dust emission, syn- 
chrotron radiation, or unresolved point sources, to men- 



tion just a few (Banday et al. 1996). One should also take 
into account cosmic variance and additional non-Gaussian 



instrumental noise in the observational data (Tegmark 
1997Q7] 

Therefore, establishing the Gaussian nature of the sig- 
nal or detecting some distinctive non-Gaussianity is cru- 
cial. In fact, it would allow us to reveal the nature of the 
primordial fluctuations. For instance, a confirmation that 
the fluctuations of the CMB are Gaussian would leave 
practically no alternative to standard inflation, since it 
would definitely rule out most of the models that predict 
non-Gaussian fluctuations and are at present still com- 
patible with observations. In addition, the investigation of 
non-Gaussianity is also important for a practical reason. 
Most of the current techniques for estimating the power 
spectrum from the observed signal with significant noise 



are optimised for the Gaussian fields only ( 


Feldman et al. 


1994 


3 


Knox et al. 199S; 


Ferreira et al. 1998 


)• 



Many authors have searched for non-Gaussian signa- 
tures in CMB data and in the large-scale structure us- 



ing such diverse tests as peak statistics (Bardeen et al 



1986: Bond and Efstathiou 1987; Vittorio and Juskicwicz 



1987j)Jthe genus curve ( poles 198$ |Smoot et al 199J ), 
higher-order correlations ( Luo and Schramm 1993| ), peak 
correlations (Kogut et al. 1996|), and global Minkowski 



functionals ( Gott III et al. 1990|; Winitzki and Kosowsky 



1997; Echmalzing and Gorski 1998) 



The first analysis of two-dimensional theoretical maps 
of the temperature fluctuations that used the total area, 
length of the boundary and genus for the excursion set 
was done by Gott III et al. (199C) although without re- 
ferring to Minkowski functionals. Schmalzing and Gorski 
(1998[]]discuss the application of the Minkowski function- 
als to all-sky maps taking into account the curvature of 
the celestial sphere. They apply these statistics to the 
high-latitude portion of the four-year COBE DMR data 
and argued for its advantage as a test of Gaussian signal. 
They conclude th at the field is cons istent with a Gaus- 
sian random field. Colley et al. (1996 ) measure the genus 
of the temperature fluctuations in the COBE DMR 4- 



year sky maps and come to a similar conclusion. Heavens 
(1999jJ] compute the bispectrum of the four-year COBE 
datasets and again fail to find evidence for non-Gaussian 
behaviour. 



However, Ferreira et al. (1998) study the distribution 
of an estimator for the normalised bispectrum and con- 
clude that Gaussianity is ruled out at the confidence level 



at least of 99%. In a recent paper, Novikov et al. (1999) 
suggest to use the partial Minkowski functionals as quan- 
titative descriptors of the geometrical properties of CMB 
maps. They apply their technique to the four-year DMR 
COBE data corrected for the Galaxy contamination and 



later shown Banday et al. (2000) that the observed non- 
Gaussianity can be explained by systematics. 

While these findings may appear confusing at first 



sight (Bromley and Tegmark 1999), they actually high- 
light the importance of studying non-Gaussian fields in 
a quantitative way. After all, while the Gaussian random 
field is a well-defined notion, non-Gaussian fields are still 



have to be substantiated. For instance, Kogut et al. (1996) 
introduce a non-Gaussian model where perturbations in t— 
space are drawn from inde pendent xt d istributions with v 
degrees of freedom, while White (1999 ) puts forward the 



use of xt distributions in real rather than in ^-space. 

In our paper, we consider a cosmological model with 
non-Gaussian initial conditions, with a statistics in real 
space following a x 2 distribution with one degree of free- 
dom. This type of perturbations results in very special 
non-Gaussianity of the CMB anisotropy. Therefore, it is 
well suited for the study of statistical methods. We calcu- 
late the Minkowski functionals and the distribution func- 
tions of extrema of the CMB field smoothed with COBE 
resolution in this model, and compare the outcome to the 
results for a Gaussian field with the same spectrum. 

This article is organised as follows. In Section |^ we 
briefly review the definition of the Minkowski functionals 
and the distribution of extrema in two dimensions, and 
compare their expectation values for a Gaussian field and 
a x 2 field with the same smoothing scale and spectrum. 
As an aside, we investigate the behaviour of the one-point 
probability distribution under smoothing. In Section |3| we 
first outline an algorithm for doing Fast Fourier Trans- 
forms on the sphere that was used for our CMB map 
simulations. Afterwards, we compare the Minkowski func- 
tionals and the peak statistics of the COBE DMR data to 
both the Gaussian and the non-Gaussian model. Section |J 
summarises and provides an outlook. 

2. Theory 

2.1. Gaussian and x 2 random fields 

Let us consider the temperature fluctuations §3- of the 
CMB temperature on the sky, parametrised with spheri- 
cal coordinates d and cp. Normalising with respect to the 

dispersion yfa, where a = ^ ( Qjh ) 2 ^> is the variance of the 

temperature fluctuations, we obtain a random field ip) 
with zero mean and unit variance: (it) = 0, (u 2 ) = 1. 

In the following, we will consider two models for the 
random field given by the normalised temperature fluctu- 
ations on the sky. 

The standard is to model u as a Gaussian random field. 
The properties of Gaussian random fields are very well 
known (see e.g. Adler 1981). 

Apart from that, we will also use a x 2 n °ld with one 
degree of freedom, as suggested by the model of Lindc and 
Mukhanov (1997). In order to retain zero mean and unit 



D. Novikov, J. Schmalzing, and V. Mukhanov: On non-Gaussianity in the Cosmic Microwave Background 3 




Fig. 1. The analytical expectation values of the 
Minkowski functionals for the Gaussian field (solid) and 
the x 2 ne ld (dashed). 



Fig. 2. The analytical expectation values of the extrema 
distributions for the Gaussian field (solid) and the \ 2 field 
(dashed). 



variance, we use the relation 



V2 



(1) 



to calculate a realisation ip of a \ 2 field from a given re- 
alisation d> of a Gaussian random field. 



2. 2. Minkowski functionals and peak statistics 

Choosing a threshold i>, we can divide the sphere into 
two parts: Hot regions where the random field u passes 
the threshold, and cold regions where u < v. The hot 
region is also called the excursion set of the field u over the 
threshold v. Its properties can be characterised by a large 
variety of geometrical quantities, including the Minkowski 



Boundary length: v\{y) is proportional to the total 
length of the boundary between cold and hot regions. 

Euler characteristic: v 2 {v), being a purely topological 
quantity, counts the number of isolated hot regions mi- 
nus the number of isolated cold regions. This is not ex- 
actly true on the sphere, where a generalised Gauss- 



Bonnet theorem holds ( Allendoerfer and Weil 1943) 



ticpa 
'ro'm- 



functio pals ( Minkowski 1902 



but nevertheless we will use this relation in the follow- 
ing calculations. 

Since Minkowski functionals are additive with respect 
to isolated parts of the sky, they can be used for patchy or 
incomplete coverage. Their calculation requires only 0(N) 
operations for calculating the Minkowski functionals on a 
map of N pixels, once the first and seco nd derivatives of 
the field at the pixel locations are known ( [Schmalzing and 



Fro m the mathematical point — of view, — Minkowski 



Gorski 1998) 



functionals have several distinguishing properties. Had 
wiger (1957) showed that all global morphological de- 
scriptors (satisfying translational invariance and additiv- 
ity) for patterns in d-dimensional space are linear combi- 
nations of just d + 1 Minkowski functionals v^(v), with 
/i G {0, ... ,d}. Moreover, in up to four dimensions all 
Minkowski functionals have simple and intuitive geomet- 
rical meanings. In two dimensions these are: 

Area: vq(is) is the total area of all hot regions, that is 
points where u(?9, if) > v. 



Both for a Gaussian and a \ random field, the 
Minkowski functionals in two dimensions arc known ana- 



lytically ( [Tomita 1990| ; |Worsley 1994j ; |Schmalzing 1999| ). 
Since the field is normalised to unit variance, the analyti- 
cal values depend on a single parameter r = (it 2 ), given as 
the variance of any of the field's first derivatives. Note that 
t has dimensions of inverse length squared, so one usually 
interprets r corr = (2t) -1 / 2 as the so-called "correlation 
length" of the random field u. 



4 D. Novikov, J. Schmalzing, and V. Mukhanov: On non-Gaussianity in the Cosmic Microwave Background 



For the Gaussian random field, we have^] 



f \ 1 J * 



2 < 




(2) 



For the \ 2 random field with one degree of freedom we 



have 



vi{v) 



4V2 
T 



I 1 Z/ \ 

exp I 1 I 



(3) 



47iyFV 2 V2 



1 v ( 1 v\ 
cxp I 1 I 



V2j 



Figure [I] shows a comparison of the expectation val- 
ues of the Minkowski functionals of the Gaussian and x 2 
random fields for the same parameter r. 

By the Morse theorem ( Morse and Cairns 1969| ), the 
Eulcr characteristic x{ u ) °f the excursion set above a cer- 
tain threshold v is related to the number of extrema above 
this threshold by 



X(v) = N max (v) + N mia {v) - N sad (v). 



(4) 



Here iV max (i/), A^ m in(^) and iV sa d(z/) denote the number 
of maxima, minima and saddle points, respectively, where 
the value of the field itself lies above the threshold v. Along 
with the Minkowski functionals, we use these three distri- 
bution functions of all possible kinds of extrema - maxima, 
minima and saddle points - as another measure of non- 
Gaussianity. We always normalise these numbers by the 
total number N ext = A r roax (-oo)+iV niill (-oo)+iV sa d(-cxD) 
of extrema in the field, giving number densities n max (v), 
n m m(v) and n sa d(^), respectively. 

The peak statistics for the Gaussian random field have 



been thoroughly investigated (Bond and Efstathiou 1987; 



Bardeen et al. 1986), but unfortunately cannot be written 



in a simple form. The same holds for the \ 2 field, whose 
extrema distribution functions are in principle easily ob- 
tained from the Gaussian case via Equation (Q). Further- 
more, the shapes of the Minkowski functional curves are 
universal, that is they do not depend on anything but 
the statistical nature (Gaussian or \ 2 ) of the underlying 
random field, and the analytical expectation values only 
contain constants of proportionality, parametrised by the 
correlation radius r corr . This is not the case for the distri- 
bution functions of extrema, where a second parameter 7 

1 $(2;) = -j= f^die - ' is the error function. 

2 The results for an arbitrary nu mber of degrees of freedom 
can be found in |Schmalzing (1998| ). 



also affects the shape of the curves. We demonstrate this 
in Figure ^|, by showing the expectation values for two 
different values of 7. 

Since the Euler characteristic is proportional to the 
third Minkowski functional v%, the distribution functions 
of extrema and the Minkowski functionals are not inde- 
pendent. In the following Section ||, we will use the first 
two Minkowski functionals vq and v\, and the three dis- 
tribution functions of extrema to test our models against 
the four-year COBE DMR data. However, let us first take 
a look at the effects of smoothing on non-Gaussian fields. 

2.3. Non-Gaussianity and Gaussian smoothing 

As an aside, we note that CMB temperature maps usu- 
ally suffer from considerable smearing through the finite 
instrument beam size and are usually smoothed in order 
to reduce the noise level. Therefore, a non-Gaussian CMB 
signal may appear in the map with a much distorted distri- 
bution. Here, we attempt to quantify this effect. We focus 
on the one-point probability distribution (PDF) of the 
field, because it is easy to handle and also closely related 
to the zeroth Minkowski functional, being nothing but its 
nagative first derivative. It is of course conceivable that 
a non-Gaussian random field still has a Gaussian PDF. 
Other statistics can in principle be treated in a similar 
fashion. 

We start from the cosmological CMB signal and con- 
sider it as a random field u(x) on some <i-dimensional 
space, where in our case d = 2. Introducing a smooth- 
ing filter gr(x, t), where t is the smoothing scale and 
g(x, 0) = 5(x) - since the unsmoothed field should be 
equal to the original field - we obtain the smoothed field 
u(x, t) by convolution: 



u{x,t)=AT(t) d d yg(x-y,t)u(y), 



(5) 



where the constant Af(t) is chosen such that the smoothed 
field remains normalised to unit variance. We may th ink 
of it(x, t) as a field in d+ 1-dimensional "scale space" ( tci 
Baar Romeny et al. 1991 ). 

If the filter g is Gaussian, it obeys the diffusion equa- 
tion 



dg(x, t) 
dt 



iA fl (x,t), 



(6) 



where the Laplacian A is of course taken with respect to 
the spatial coordinates x only. 

Combining Equations (^|) and (|^), we can write down 
an "evolution equation" for the field u as the scale t 
changes: 



<9w(x, t) 
dt 



t (A 



u(x, t). 



(7) 



The second expression r~^ VI in the operator enters through 
the scale-dependence of the normalisation factor JV(t). 
Details of the calculation can be found in Appendix |a|. 
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Fig. 3. The PDFs of a smoothed x 2 random field (solid lines). From left to right, the smoothing scale increases, 
ranging from no smoothing (left panel) over little smoothing (middle panel) to considerable smoothing (right panel). 
The transition of the curve's shape from clearly non-Gaussian to almost Gaussian is obvious. For comparison, the 
Gaussian bell curve is plotted in all three examples (dashed lines). 



This equation enables us to study the one-point prob- 
ability distribution P(u,t) of the smoothed field u(-,t) as 
the smoothing scale changes. Writing this probability den- 
sity as 

P(u,t) = (5(u(x,t) - u)) (8) 

and applying the partial derivative with respect to t, we 
immediately obtain 



dP(u,t) 



8 



(Au) u 



P(u,t) 



(9) 



dt du 

We refer the reader to Appendix |B| for the technical de- 
tails of this calculation. The quantity (Au) u denotes the 
average of the Laplacian of the field it(x, t) under the con- 
dition that its value u is fixed. It is interesting to note 
that this equation is written in conservative form, that is 
its integral overdu vanishes. 

For a variety of random fields, the conditional aver- 
age of the Laplacian can be calculated analytically. Most 
notably, for both the Gaussian and the \ 2 field (Am)„ = 
— . In these cases, the right-hand side of Equation (^|) 
simply becomes equal to zero. 

As far as the Gaussian random field is concerned, 
P(u) cx exp(— u 2 /2) is a stationary solution of Equa- 
tion (^|) . This reflects the well-known fact that a Gaussian 
random field stays Gaussian under smoothing, and in fact 
under any linear filtering. 

In the case of the \ 2 field, however, the probability dis- 
tribution is non-zero only for u < 1/V2 by Equation ([!]), 
and is not differentiable at the upper bound of the field. 
Therefore, under smoothing the field evolves away from 
the x 2 distribution, and finally ends up with a distribution 
close to the stationary Gaussian solution at large smooth- 
ing lengths. Figure || illustrates the effect using a x 2 ran- 
dom field with various degrees of smoothing. 



This shows that although our two models are fairly dif- 
ferent at first glance, the noise, beam size and other effects 
of actual observations tend to make the fields appear more 
similar. We shall look upon the practical consequences of 
this now, using the four-year COBE DMR data as an ex- 
ample. 

3. Data analysis 

We illustrate our methods on the COBE DMR four-year 
all-sky maps ( Bennett ct al. 1 99(f ). We work with the map 
constructed from all three bands and both channels of the 
DMR instrument by the so-called subtraction method. 
Although this method is supposed to remove most of the 
galactic contamination, we avoid the galactic plane com- 
pletely by applying a cut to 45° latitude. 

For the statistical tests presented in the following, we 
use a large number of realizations of the Gaussian and 
X 2 random fields. In both cases, we calculate the fields 
using a scale-free spectrum with a correlation radius of 
1°. To mimic the DMR beam, we smooth the fields with 
a Gaussian filter of 7° FWHM. Examples of realisations 
of the Gaussian and the x 2 random fields, together with 
the COBE DMR maps, are shown in Figure 0. 



3.1. Fast Fourier Transform on the sphere 

Statistical tests on CMB maps usually involve simulat- 
ing large numbers of model maps from given power spec- 
tra. The naive approach to this computational task is ex- 
tremely time-consuming, requiring 0(N 2 M 2 ) operations, 
where N and M are the number of pixels on the equator 
and the number of multipoles, respectively. Here we use a 
Fast Fourier Transform (FFT) in flat space over spherical 
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COBE north pole COBE south pole <p [degrees] 



Fig. 4. Four maps showing realisations of the Gaussian 
and x 2 random fields and the northern and southern parts 
of the COBE DMR map, respectively. All maps extend 
from the pole to 30° latitude. The lines show isocontours 
of the field at half integer multiples of the standard devi- 
ation. Shaded areas indicate regions where the field value 
drops below the mean. Open circles show the positions 
of minima, full circles correspond to maxima, and saddle 
points are indicated by drude's feet. 



coordinates outlined by Risbo (1996| ). This means that the 
FFT is used in both the polar and the azimuthal direction, 
while other methods usually perform the FFT over the po- 
lar angle alone. Let us first outline the mathematical and 
computational ideas. 

Given a realisation of a random field in terms of its 
spherical harmonics coefficients ae m , we wish to calculate 
its values u(#, <p) in real space parametrised with spherical 
coordinates (z?, ip). This can be done through the spherical 
harmonics expansion 



m-i +e 

E E 



at m Ytm(&,<p), 



(10) 



where the functions Yi m are the spherical harmonics. They 
are related to the associated Lcgcndre functions Pi m by 



Y lm (d, v ) = h lm P em (cos$)e mip , 



(11) 



where m and £ are called the order and the degree of the 
function, respectively, and 



>2l+l {i-m)\ 
4tt (£ + m)\ 



(12) 



Fig. 5. A COBE all-sky map, projected onto the plane 
in cylindrical coordinates and mirrored at the North pole 
to give a square suitable for the FFT method described 
in the text. The meaning of the various line and symbol 
styles is explained in the caption of Figure 0. 



All spherical harmonics Yg, m enjoy the symmetry prop- 



erty 



Y tm (0,<p) = Y lm (2Tr-d,n + ip). 



(13) 



This can be used to extend the field ip), which is nor- 
mally defined for (??, ip) E [0, 7r] x [0, 2ir], to the whole 
square [0, 2ir] 2 , simply by mirroring the field at the point 
(n, 7t). The resulting square map is illustrated by an ex- 
ample in Figure |^. 

The extended field in real space, which now consists 
of two identical copies of the actual map, is periodic both 
in $ and in ip, and can therefore be expanded into a real- 
space Fourier series: 



M-l M-l 

i V>*>)= E E bjkJ'****, 

j=0 k=0 



(14) 



Inverting Equation ( pL4|) and inserting Equation (10) 
allows us to establish a relation between the Fourier coef- 
ficients and the spherical harmonics coefficients ai m '- 



M-l 



bjk — J]] aimhl m Phnj\ m=k 
M-l 

^2 a f^ h imPlraj\ m=k _ M , (15) 



£=M-k 
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P 

area 


-^length 


p 

i max 


1 mm 




South 


Gaussian 


94% 


12% 


38% 


52% 


53% 




xl 


87% 


61% 


44% 


56% 


57% 


North 


Gaussian 


33% 


44% 


89% 


87% 


95% 




xl 


31% 


28% 


76% 


81% 


88% 



Table 1. Summary of the probabilities of obtaining the 
COBE DMR data as a realisation of the models - a Gaus- 
sian random field and a \ 2 random field - considered 
here. The analysis is done separately for the Northern and 
Southern polar caps. Since we use five related statistics to 
obtain these probabilities, five values are shown for each 
case. 



where the pe m j are the coefficients of the Fourier series of 
the associated Legendre functions: 



£m 



cost?) 



M-l 

E 

3=0 



e J pemj- 



(16) 



Using the well-known (Press et al. 1987) recurrence 
relation for the associated Legendre functions, it is easy 
to derive recurrence relations for numerical evaluation of 
the p tmj . 

Let us now summarise the steps for efficiently calculat- 
ing a CMB map in real space from its spherical harmonics 
coefficients. 

1. Obtain a realisation of the random field in spherical 
harmonics ag m . 

2. Calculate the intermediate Fourier coefficients bjk us- 
ing Equation (p"oj). This step requires 0(M 3 ) opera- 
tions. 

3. Perform the two-dimensional Fast Fourier Transform 
from Equation ( flil) in order to obtain the real-space 
realisation u(i), ip). This step requires 0(N 2 log M) op- 
erations. 

Obviously, our approach produces a cylindrical pixeli- 
sation in real space, resulting in non-uniform pixels. This 
makes it unsuitable for applications that involve quadra- 



tures, where uniform pixelisations (e.g. Gorski et al. 1999) 
are desirable. However, the method is well suited for the 
calculations presented below, because it allows large over- 
sampling of the map in real space, that is TV -C M, with 
little computational effort. This is desirable both for ac- 
curate calculations of Minkowski functionals of a map by 
tracing the contours with high accuracy, and for determin- 
ing the exact positions of extrcma. 

3. 2. Minkowski functionals 

The Minkowski functionals of the COBE DMR maps are 
displayed in Figure ^. Since the Northern and Southern 
polar caps were analysed separately, we show two curves 
in each panel. The shaded areas indicate the average and 



v„ - area 



Vj - length 



0.8 
0.6 
0.4 
0.2 





E u le r e h a r a e I . e r i y I . i c 



0.1 




Fig. 6. The Minkowski functionals of the COBE DMR 
data in comparison with the average values for a Gaussian 
random field. The areas indicate the standard deviation. 

standard deviation of 1000 realisations of a Gaussian ran- 
dom field with the same two-point characteristics as the 
COBE data. 

In order to assess the probability that the COBE DMR 
maps are indeed a realisation of our toy models, we use a 
non-parametric test (Novikov et al. 2000). 

We computed N = 1000 realisation of each random 
field. The nth realisation yields Minkowski functionals 
V^ n \u), with the threshold v ranging from —5 to 5. We 
calculate the average Minkowski functionals 



1 



N 



71=1 



(17) 



for the random field, and determine the Li-distance of 
each realisation from this average: 



AW = / du 



(18) 



We then obtain the same quantity A^ OBE from the data, 
and test the hypothesis that this A^ OBE is distributed 

according to the distribution function of the A^ of the 
model we are currently investigating. 

The resulting probabilities are displayed in Table [1} 
Note that we only show the values for the first two func- 
tionals, the area Vq and the length Vi, since by Equa- 
tion (0) the last Minkowski functional is a linear combi- 



nation of the extrema statistics analysed in Section 3.3 
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Fig. 7. The extrema statistics of the COBE DMR data in 
comparison with the average values for a Gaussian random 
field. The areas indicate the standard deviation. 

3. 3. Distribution of extrema 

For the COBE DMR data, the distribution functions of 
all three kinds of extrema are shown in Figure ^. For com- 
parison, we also show the mean and standard deviations 
over 1000 realisations of a Gaussian random field with the 
same two-point characteristics. 

In order to quantify the difference of the COBE data 
and t he models, we use a standard Kolmogorov-Smirnov 
test ( Press et al. 1987 ). The resulting probabilities are 



shown in Table |lj. 
4. Summary and Outlook 

We have investigated a method for probing the non- 
Gaussianity of CMB maps, based on Minkowski function- 
als and on distribution functions of extrema. In order to 
make the deviations from Gaussianity quantitative, we 



used a model proposed by Linde and Mukhanov (1997), 
producing to a % 2 distribution with one degree of freedom 
for the primordial fluctuations. Analytical calculations for 
our statistics showed that the non-Gaussian and Gaussian 
models produce markedly different results. 

The situation may change, however, when smoothing 
is taken into account. We derived an equation govern- 
ing the change of the one-point probability distribution 
with the smoothing scale and found that, while the Gaus- 
sian field remains Gaussian under smoothing, the x 2 field 
evolves towards a more symmetric distribution with closer 



resemblance to the Gaussian bell. The outcome of this 
effect became clear in an application of the method to 
the COBE DMR data - we saw a slight preference of the 
Gaussian model, but this was by no means significant. 
Since all realisations of the maps were highly oversam- 
pled in real space, the failure to discriminate between the 
two models is definitely due to the rather large smoothing 
scale. Therefore, an application of the method to the re- 



cently complete d Boomerang flde Bernadis et al. 2000D 
and MAXIMA ( Hanany et al. 200C ) experiments appears 
highly promising, as soon as the data become available. 
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Appendix A: Derivation of Equation (|7j) 

We start from the smoothing prescription in Equation (|J) : 



(x,t) = Af(t) d d y 9 (x-y,() H (y), 



(A.l) 



and apply the partial derivative with respect to the 
smoothing scale t. Straightforward differentiation yields 



du (x, t) 
dt 



dt 



& d yg{x-y,t)u{y) + 



/d' / g(x - y '*W (A.2) 



We use the diffusion equation (g) on the second addend, 
and eliminate the integrals by reinserting the smoothing 
prescription (§). This leads to: 



dt Af(t) 



u(x,t) +tAu(x,t). (A.3) 



In order to evaluate the partial derivative of the normali- 
sation factor Af(t) , we use the normalisation condition 



(A.4) 



and again apply the partial derivative with respect to t. 
Together with Equation (|A.3[), this becomes 



= J d d xu(x,t) 



du(x,t) 
dt 



= jfhv J d d xu{x,t) 2 +t J d d xu(x,t)Au(x,t) (A.5) 

dAf(t) 
— — — t Y r~ 2 

M{t) corr ' 

The first integral is eq ual to one because of the normal- 
isation condition (A.4), while the second integral can be 
rewritten as the variance of the gradient of u{x, t) using 
partial integration, which is equal to the inverse square of 
the correlation radius r corr by definition. 

Now we can eliminate the explicit reference to the nor- 
malisation factor JV(t) from Equation (A.3) and obtain 
our final result, Equation (|t]): 



du(x,t) 
dt 



-H A + r c"orr «(X,*). 



(A.6) 



It is a common bad habit among physicists to use the 
same notation for utterly different things. We are no ex- 
pection to this, employing the same letter both for the 
smoothed random field and for the argument of the proba- 
bility distribution. However, we always add the scale space 
coordinate (x, t) to the random field u(x,t) , while the 
value u has no such argument. 

At first, the partial derivative of this equation yields 
a functional derivative of the 5-function with respect to 
the field it(x, t). Fortunately, it can be rewritten as an 
ordinary partial derivative with respect to the value u. 
Furthermore, the value u (not the field ii(x, t)\) only enters 
through the J-function. We can then insert Equation ^ 
for the evolution of the field under smoothing. All in all, 
the calculation proceeds as follows: 



dP{u,t) 
dt 



d6(u(x, t) — u) du(x,t) 

x ~bT~ 

d I du(x,t) 
--{5(u(x,t)-u)—^- 

-t— {8(u{x, t) - u) (A + r c - 2 rr ) u(x, t)) 

-*7r ( *) - u ) ( Au ( x ' *) + -r 



-t 



du 
d_ 
du 



(Au) u 



P(u,t) 



(B.2) 



so we end up with Equation 



Appendix B: Derivation of Equation (^) 

This time, we investigate the properties of the one-point 
probability distribution function P(u) under smoothing. 
We start from the definition (B): 



P(u,t) = {S(u{x,t)-u)). (B.l) 



